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■ We present RIz photometry of four consecutive transits of the newly discov- 
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i-H ■ to be .Rg = 0.928^o:ol3 #o, assuming a stellar mass of M s = 1.00 ± 0.03 M . 

^ ■ The uncertainties in the planetary and stellar radii are dominated by the uncer- 

tainty in the stellar mass. These uncertainties increase by a factor of 2-3 if a 
t^. ■ more conservative uncertainty of 0.10 M & is assumed for the stellar mass. Our 
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1. 



Introduction 



An exoplanetary transit is a rare opportunity to learn a great deal about both the 
planet and the star. With precise measurements of the amount of light blocked by the 
planet as a function of time, it is possible to infer the relative sizes of the star and planet, 
the orbital inclination, and the stellar limb-darkening function. Coupled with measurements 
of the time-variable Doppler shift of the star and an estimate of the stellar mass, one learns 
the planetary mass and the stellar radius. These fundamental measurements set the stage 
for a host of more subtle measurements of effects such as planetary atmospheric absorption 
lines, thermal emission, spin-orbit alignment, and timing anomalies, as reviewed recently by 
Charbonneau et al. (2006a). 

For these reasons, newly-discovered transiting exoplanets are welcomed with open arms. 
The tenth such object was recently reported by McCullough et al. (2006). The parent 
star, XO-1, is bright (V = 11, Gl V), making it a favorable target for precise observa- 
tions. The planet has an orbital period of ~4 days and a mass and radius comparable to 
Jupiter, although McCullough et al. (2006) point out that their photometry actually implies 
a mean density that is somewhat smaller than theoretical expectations for "hot Jupiters." 
If confirmed, this would put XO-lb in the same category as the anomalously large planet 
HD 209458b, and may have implications for the various theories that have been espoused 
for that object. 

Two of us (M.J.H. and J.N.W.) have initiated the Transit Light Curve (TLC) Project, 
a long-term campaign to build a library of high-precision transit photometry, with the dual 
goals of (1) refining the estimates of the physical and orbital parameters of the target systems, 
and (2) searching for secular and short-term variations in the transit times (and light curves) 
that would be indicative of perturbations from additional bodies (Agol et al. 2005; Holman 
& Murray 2005). Here, we present results for XO-lb that were obtained as part of this 
program. We describe the observations and the data reduction procedures in § 2. In § 3 we 
describe the model and techniques we used to estimate the physical and orbital parameters 
of the XO-1 system, and in § 4 we summarize our results. 



We observed four consecutive transits of XO-lb. According to the ephemeris provided 
by McCullough et al. (2006), 



2. 



The Observations and Data Reduction 



T C (E) = 2,453,808.9170 [HJD] + E x (3.941534 days), 



(1) 
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these transits correspond to epochs 17 through 20. We employed three different telescopes: 
the FLWO 1.2m telescope (for E = 19,20); the Palomar 1.5m telescope (for E = 17), and 
the TopHAT 0.26m telescope (for E = 17, 18, 19). 

2.1. FLWO 1.2m z Photometry 

We observed the E = 19, 20 transits (UT 2006 May 28 and June 1) with KeplerCam on 
the 1.2m (48 inch) telescope of the Fred L. Whipple Observatory (FLWO) on Mt. Hopkins, 
Arizona. This camera (P.I. D. Latham) was built for a photometric survey of the target field 
of the Kepler satellite mission (Borucki et al. 2003). It has a single 4K x 4K Fairchild 486 
CCD with a 23'1 x 23'.1 field of view. We used 2x2 binning, for which the readout and 
reset time is 11.5 s and the typical read noise is 7 e~ per binned pixel. The response of each 
amplifier deviates from linearity by less that 0.5% over the range of counts from the faintest 
to brightest comparison star. We observed through the SDSS z filter, the reddest available 
band, in order to minimize the effects of limb-darkening and color-dependent atmospheric 
extinction. The effective bandpass at the red end was limited by the quantum efficiency of 
the CCD, which drops from ~ 100% at 7500 A to ~ 10% at 10500 A. We defocused the 
telescope slightly in order to enhance the duty cycle and average over pixel-to-pixel sensitivity 
variations. The full-width at half-maximum (FWHM) of a stellar image was typically ~3 
binned pixels (2"). We used automatic guiding to maintain the locations of XO-1 and its 
comparison stars to within a few pixels over the course of both nights. On each night, we 
repeatedly took 30 second exposures for approximately 5 hr bracketing the predicted transit 
midpoint. The conditions on UT 2006 May 28 appeared photometric, and the images were 
taken through airmasses ranging from 1.00 to 1.22. The conditions on UT 2006 June 1 were 
nearly photometric, except for very thin, high clouds that passed through the field between 
UT 05:50 and UT 06:30. The airmass range on this night was from 1.00 to 1.31. 

The images were calibrated using standard IRAF 1 procedures for the overscan correc- 
tion, trimming, bias subtraction, and flat-field division. We did not attempt to correct the 
fringing that was apparent with the z filter. The fringing had a small amplitude and little 
effect on the final photometry, given the accuracy of the automatic guiding. We excluded 
three images the first night that showed significant large-scale features that were not cor- 
rected by the flat-field (presumably from clouds), and one from the second night that was 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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contaminated by an artificial satellite trail. We then performed aperture photometry of XO- 
1 and 4 nearby stars of comparable brightness and color (stars 1, 3, 6, and 7 from Table 1 
of McCullough et al. 2006), using an aperture radius of 6.5 pixels (4'.'3) for both nights. We 
subtracted the underlying contribution from the sky, after estimating its brightness within 
an annulus ranging from 30 to 35 pixels in radius, centered on each star. We divided the 
flux of XO-1 by the total flux of the comparison stars. We then fit a linear function of time 
to the pre-ingress and post-egress data, and divided the entire time series by this function, 
in order to normalize the flux and remove residual systematic effects. A function of time 
proved to be a slightly better fit than the more traditional function of airmass. 

To estimate the uncertainties in our photometry, we computed the quadrature sum of 
the errors due to Poisson noise of the stars (both XO-1 and the comparison stars), Poisson 
noise of the sky background, readout noise, and scintillation noise (as estimated according 
to the empirical formulas of Young 1967 and Dravins et al. 1998). The dominant term is 
the Poisson noise from XO-1. The final time series is plotted in Fig. 1 and is available in 
electronic form in Table 1. 
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Fig. 1. — Relative z band photometry of XO-1. The left panels show the photometry (points) 
and the best-fitting model (solid line). The right panels show the residuals (observed — cal- 
culated) and representative 1 a error bars that have been rescaled such that the \ 2 P er degree 
of freedom for each time series is unity. For the first transit (E = 19), the RMS residual is 
0.15%. For the second transit, the RMS residual is 0.12%. 
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2.2. TopHAT I Photometry 

We used TopHAT to observe the E = 17, 18, 19 transits of XO-lb (UT 2006 May 20, 
24, and 28). TopHAT is an automated telescope, also located on Mt. Hopkins, Arizona, 
which was designed for photometric follow-up observations of transiting exoplanet candidates 
identified by the HAT Network (Bakos et al. 2004). It consists of a 0.26m diameter f/5 
commercially-available Baker Ritchey- Chretien telescope on an equatorial fork mount. A 
l?25-square field of view is imaged onto a 2Kx2K Peltier-cooled, thinned CCD detector, 
yielding a pixel scale of 272. In order to extend the integration times and increase the duty 
cycle of the observations, we applied a slight defocusing. The resulting PSF had a FWHM 
of 2.1 pixels (4'/6). On each night, we observed for approximately 5 hrs. 

We calibrated the images by subtracting the overscan bias and a scaled dark image, and 
dividing by an average sky flat from which large outliers had been rejected. We performed 
aperture photometry on XO-1 and on an additional ~800 stars in the field, using an aperture 
of radius 5 pixels (ll'.'O), and an exterior annulus for sky subtraction. Most of the 800 stars 
(after removing variables) were used as calibrators in a statistically-weighted manner to 
transform the magnitudes of the individual frames to the instrumental magnitude system 
of a selected reference frame. The derived light-curve still suffers from small-amplitude 
systematic errors. In order to minimize them, we used all the out-of-transit data (assuming 
it to be of constant brightness) to find the correlation with the airmass, hour-angle (linear 
fits), the pixel position (sinusoidal function), and the Gaussian profile parameters (second 
order fits). The fitted function was then applied to and subtracted from the entire light- 
curve, including the transits. These corrections in this post-processing step were of the order 
of 3 mmag or less for unsaturated points. The resulting time series is shown in Fig. 2 and 
listed in Table 1, along with the Palomar data described below. 
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Fig. 2. — Relative R and / band photometry of XO-1. The left panels show the photometry 
(points) and the best-fitting model (solid line). The right panels show the residuals (ob- 
served — calculated) and representative 1 a error bars, estimated as described in the text. 
From top to bottom, the RMS residuals are 0.21%, 0.28%, 0.31%, and 0.32%. 
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2.3. Palomar 1.5m R Photometry 

We observed the E — 17 transit with the 1.5m (60 inch) telescope at Palomar Obser- 
vatory. The CCD camera has 2Kx2K pixels with a plate scale of 0'/378 pixel -1 . In order 
to increase the duty cycle, we read out only half of the available field of view, yielding an 
effective field size of 12!9 x 6'5. The sky was cloud-free, the typical seeing was 1'.'5, and the 
observations ranged in airmass from 1.01 — 1.45. We gathered 209 i?-band images spanning 
5.2 h, with integration times of 20 s and a cadence that increased from 65 s to 90 s over 
the course of the observing sequence. The observing sequence was briefly interrupted twice 
(once prior to ingress, and once near mid-transit) by telescope calibration scripts that are 
required as part of its robotic operation. 

We used the automated P60 reduction pipeline to calibrate the images. This pipeline 
trims the overscan columns, subtracts the bias level, divides by a dome flat, and flags bad 
pixels. We performed aperture photometry of XO-1 and 5 comparison stars, using an aper- 
ture of radius 9 pixels (3'/4), and an annulus for sky subtraction ranging in radius from 30 
to 45 pixels. We divided the flux of XO-1 by the sum of fluxes of the comparison stars, 
and normalized the resulting time series to produce a flux of unity in the pre-ingress data. 
The RMS variation of the pre-ingress data is 0.17%, which we adopted as the photometric 
uncertainty in each data point. The pointing drifted by as much as 35 pixels (13") during 
the course of the observations. We believe that this drift, coupled with uncertainties in the 
flat-field image, is responsible for the poorer quality of the Palomar 1.5m data compared to 
those from the FLWO 1.2m. The Palomar data are presented in Fig. 2 and listed in Table 1. 

3. The Model 

The planetary, stellar, and orbital parameters were inferred by fitting a parameterized 
model to all of the photometry simultaneously. The model is based on a star and a planet 
on a circular orbit about their center of mass. 2 The star has a mass Mg and radius R$, and 
the planet has a mass Mp and radius Rp. The orbit has a period P and an inclination i 
relative to the sky plane. We define the coordinate system such that 0° < % < 90°. The 
initial condition is specified by T c , a particular time of conjunction (the transit midpoint). 
When the planet is projected in front of the star, the model flux declines by an amount 
that depends on the projected separation, on the stellar limb-darkening function, and also 



2 We assume a circular orbit because the radial velocity data show no evidence for a non-zero eccentricity, 
and in the absence of any evidence for additional bodies in the system, it is expected that tidal interactions 
have had sufficient time to circularize the orbit (see, e.g., Rasio et al. 1996; Trilling 2000; Dobbs-Dixon et al. 
2004). 
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on the planet-to-star area ratio. To compute this flux decrement, we assume a quadratic 
limb-darkening law, and employ the analytic formulas of Mandel & Agol (2002). The R 
and / band data are not of sufficiently high signal-to- noise ratio to justify this treatment; 
instead, we fix both combinations of the R and / limb-darkening parameters at the values 
estimated by Claret (2004) for a star of the observed effective temperature, surface gravity, 
and metallicity (see Table 5 of McCullough et al. 2006). 

We allow the T c for each of the 4 transits to be an independent parameter. This is 
because we seek to measure or bound any timing anomalies that may indicate the presence 
of moons or additional planets in the system. Obviously if we allow each of the four T c 
values to vary, we cannot independently determine the orbital period. Instead, we fix P = 
3.941534 days, the value reported by McCullough et al. (2006). This mean period is based 
on observations spanning a few years and is known more accurately than we could hope to 
determine from our time baseline of 12 days. The quoted uncertainty in the mean period is 
only 0.000027 days and is negligible for our purposes. 

There is a well-known degeneracy among M s , R s and R P that prevents all three pa- 
rameters from being uniquely determined from transit photometry alone, unless a stellar 
mass-radius relation is assumed (Seager & Mallen-Ornelas 2003). We fix M s = 1.0 M©, 
based on the spectroscopic estimate by McCullough et al. (2006). Our results may be scaled 
to other choices for the stellar mass according to R s oc (M s /M ) 1/3 and R P oc (Ms/Mo) 1 / 3 . 
The planetary mass Mp is irrelevant to the model except for its minuscule effect on the 
relation between P and the semimajor axis; for completeness, we assume Mp = 0.9 Mj up , 
again following McCullough et al. (2006). 

In total, there are 8 free parameters describing 1309 photometric data points. The free 
parameters are Rs, Rp, and i; the z-band limb darkening parameters Ui and Ui] and the 
four values of T c . In practice, we found it better to fit for the parameters 2u\ + and 
u\ — 2u2, because the resulting uncertainties in those parameters are uncorrelated (as will 
be shown below). We allowed the limb darkening parameters to range only over the values 
that produce a monotonically decreasing intensity from the center of the star to the limb. 

Prior to fitting the full set of observations, we fitted each of the six time series separately 
and determined the minimum x 2 in each case. The resulting values of x 2 P er degree of 
freedom were: 2.03 (FLWO 1.2m night 1), 1.39 (FLWO 1.2m night 2), 1.53 (Palomar 1.5m), 
0.25 (TopHAT night 1), 0.50 (TopHAT night 2), and 0.34 (TopHAT night 3). Thus, it seems 
that the calculated uncertainties were somewhat underestimated for the FLWO and Palomar 
data, and overestimated for the TopHAT data. Before proceeding, we scaled the estimated 
uncertainties of each time series individually so that the resulting value of x 2 P er degree of 
freedom was unity. Table 2 gives the uncertainties after this scaling was performed. 
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We determined the best-fitting model using the AMOEBA algorithm (Press et al. 1992) 
to minimize the error statistic 



where /j(obs) is the flux observed at time tj, <jj is the corresponding uncertainty, and /j(calc) 
is the calculated value. In Figs. 1 and 2, the left-hand panels show the best-fitting model as 
a solid line, and the right-hand panels show the results of subtracting the calculated values 
from the observed values. The FLWO z band data show random-looking residuals with a 
standard deviation of 0.15% for the first night, and 0.12% for the second night. Almost all of 
the leverage on the stellar and planetary parameters comes from these data. The P60 data 
also show nearly random residuals, but at a higher level of 0.2%, and with occasional outliers. 
The TopHAT data are noisier, with a standard deviation of 0.3%, and show some signs of 
correlated residuals (i.e. systematic errors). The uncertainties in the fitted parameters were 
estimated using two different methods, described below. 

The first method was a bootstrap analysis, similar to those we have performed previously 
for the transiting exoplanets HD 209458b (Winn et al. 2005), OGLE-TR-10 (Holman et al. 
2005), and HD 149026b (Charbonneau et al. 2006b). We refitted the parameters to each of 
10 4 different "realizations" of the data. These realizations were sets of 1309 data points drawn 
randomly from the actual data set, with duplications allowed (i.e., with replacement). Each 
realization was required to preserve the total number of points in each of the 6 individual 
time series. The resulting collection of 10 4 optimized parameter sets was taken to be the 
joint probability distribution for the parameters. 

The second method was a Markov Chain Monte Carlo (MCMC) simulation. In this 
method (described lucidly for astrophysicists by Tegmark et al. 2004 and Ford 2005), a 
stochastic process is used to create a sequence of points in parameter space that approxi- 
mates the desired probability distribution. The sequence, or "chain," is generated from an 
initial point by iterating a "jump function." In our case the jump function was the addition 
of a Gaussian random number to each parameter value. If the new point has a lower x 2 than 
the previous point, the jump is executed; if not, the jump is only executed with probabil- 
ity exp(— A% 2 /2). Under fairly benign mathematical assumptions, the chain will eventually 
converge to the desired probability distribution. To speed convergence, the Gaussian per- 
turbations should be large, but not so large that all jumps are rejected. We set the relative 
sizes of the perturbations using the 1 a uncertainties estimated previously by the bootstrap 
method, and we set the overall jump size by requiring that ~25% of jumps are executed. 
We created 10 independent chains, each from a random initial position ~5o" away from the 
optimized parameter values. Each chain had 500,000 points, the first 20% of which were dis- 





- 11 - 



carded to minimize the effect of the initial condition. The typical correlation length for each 
parameter (see Tegmark et al. 2004) was ~400 points, giving an effective length of ~1000 
per chain. To check the convergence and the consistency between chains, we computed the 
Gelman & Rubin (1992) statistic for each parameter, which is a comparison between the 
inter-chain variance and the intra-chain variance. The results were within a few per cent of 
unity, a sign of good mixing and convergence. 

The two methods produced very similar results. We also checked the results for the 
planetary and stellar radii using the more traditional approach of stepping each parameter 
through a sequence of values while allowing all of the other parameters to float, and iden- 
tifying the values for which Ax 2 = 1 as the 68% confidence limits (as done by Brown et 
al. 2001, among others). Again, the values and the uncertainties were comparable to the 
results of the MCMC and bootstrap methods. For brevity, we report only the MCMC results 
for the remainder of this paper. The probability distributions for some of the parameters are 
shown in Fig. 3, and some of the correlations between the parameters are shown in Fig. 4. 
Table 2 lists the median value of each parameter with their 68% confidence limits, based on 
the MCMC results. 
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Fig. 3. — Estimated probability distributions of some planetary, stellar, and orbital param- 
eters. The histograms show the results of 10 Markov Chain Monte Carlo simulations, each 
with 400,000 points. The median of each distribution is indicated with a solid line. The 
dashed lines enclose 68% of the results, with equal probability on either side of the the me- 
dian. The arrows show the the choice of parameters that minimizes x 2 - The numbers in 
Table 2 are the median values, with confidence limits given by the dashed lines. 
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Fig. 4. — Joint probability distributions of some planetary, stellar, orbital, and limb- 
darkening parameters, based on the MCMC simulations. The contours are isoprobability 
contours enclosing 68% and 95% of the points in the Markov chains. In the Ui-u 2 plot, the 
dot shows the values calculated by Claret (2004). In the other plots involving limb darkening, 
the dotted lines shows the values calculated by Claret (2004). 
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4. Summary and Discussion 

Through observations of four consecutive transits, we have significantly improved upon 
the estimates of the system parameters of XO-1. The most interesting parameters are the 
radius of the star, the radius of the planet, and the mid-transit times, which will be discussed 
shortly. The results for the other parameters are not especially interesting but they do seem 
reasonable. The results for the orbital inclination are best described as bounds on the 
impact parameter b, which is the minimum projected star-planet distance, in units of the 
stellar radius. It is given by b = acosi/Rs, where a is the orbital distance. The data favor 
a central transit, with b < 0.27 and % < 88?53 at the 95% confidence level. 

Although the survey and follow-up photometry of McCullough et al. (2006) were im- 
pressive, and built a convincing case for an exoplanet, those authors did not attempt to fit 
for the stellar radius when modeling the transit light curve. Instead, they used the value 
Rs/Rq = 1.00 ± 0.08, based on an interpretation of the stellar spectrum. This is because 
the inference of R$ from a transit light curve requires that the ingress and egress are well 
sampled and measured with a high signal-to-noise ratio. This type of data was not available. 
The higher precision and finer time sampling of our data, and of the z band data in par- 
ticular, allow for the determination of Rs from the light curve, without relying on spectral 
modeling and theoretical isochrones. The resulting "photometric" value of R s is still subject 
to a systematic error due to the covariance with M s , but the dependence is fairly weak, 

1/3 

Rs oc M s , generally leading to a smaller uncertainty in Rs than can be achieved from 
spectral modeling and theoretical isochrones. 

Our result is Rs/R G = 0.928t2;2]J, which is consistent with (but more precise than) 
the value determined by McCullough et al. (2006). Here we have incorporated the 0.03 M Q 
uncertainty in Ms determined by McCullough et al. (2006). We note that this radius is 
somewhat small for the Gl V spectral type of XO-1, but it is still consistent, given the stated 
uncertainties. We remind the reader again that the quoted result assumes Ms = 1.0 M Q , 
and that the inferred Rs scales as (Ms/M Q ) 1//3 . From the Yonsei-Yale isochrones, a stellar 
mass of M s = 0.96 M corresponds to a radius of Rs = 0.91 R Q , for solar metallicity and an 
arbitrary age of 3.6 Gyr, (Yi et al. 2001). Thus, a 1.3 a change in the estimated stellar mass 
yields an estimated stellar radius that is precisely in line with theoretical expectations. We 
also note that the stellar radius uncertainty is a factor of 2-3 larger if a more conservative 
uncertainty of 0.10 M Q is assumed for the stellar mass, as shown in Table 2. 

Our derived radius of XO-lb is R P /Rj U p = 1.184±g;g?| (again assuming the uncertainty 
in the stellar mass to be 0.03 M ). Previously, McCullough et al. (2006) found Rp/Rj up = 
1.30 ±0.11. These figures are also in agreement right within their respective 68% confidence 
limits. Interestingly, we obtain very precise agreement with McCullough et al. (2006) for all 
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parameters if we first time-average our data into 5-minute bins (i.e., by a factor of 8, for 
the z band data). The McCullough et al. (2006) data were averaged into bins ranging in 
width from 3 to 9 minutes, depending on the telescope used. We suggest that it is possible 
that some of the previous results were slightly biased by the coarser time sampling of the 
photometry. Resolving the degeneracy among the stellar radius, planetary radius, and orbital 
inclination requires adequate sampling of ingress and egress. 

The uncertainties in the limb darkening parameters U\ and u 2 are highly correlated, 
with the linear combination 1u\ + u 2 being well constrained by the data, and the orthogonal 
combination u\ — 2u 2 being weakly constrained by the data (see the lower left panel of Fig. 3). 
We find 2ui + u 2 = 0.86 ± 0.05. This is 2 o larger than the value based on the theoretical 
calculations of Claret (2004), which predict u x = 0.21, u 2 = 0.33, and 2u x + u 2 = 0.75 (for 
the standard z band, T = 5750 K, log g = 4.5, [M/H\ = 0.05, and microturbulent velocity 
v t = 2.0 km/s). The theoretical values are shown in Fig. 3 as the solid symbol in the u\-u 2 
plot, and as dotted lines in the other two limb-darkening plots. One might consider reducing 
the number of degrees of freedom and adopting the Claret (2004) values as fixed quantities. 
When we do so, we find R S /R Q = 0.94, R P /Rj up = 1.22, and b = 0.26 (i = 88.65 deg), 
with the minimum x 2 increased by 9. However, given the quality of the z band data, the 
unknown level of uncertainty in the theoretical values, and the possible differences between 
the FLW048/KeplerCam z band and the standard SDSS z band, we believe fitting for the 
limb-darkening coefficients is more appropriate. 

Some of the probability distributions shown in Fig. 3 are asymmetric. This is typical 
of all fits to transit light curve data. The fact that the orbital inclination has a maximum 
value (namely, 90°), combined with the measured durations of the ingress, egress, and the 
full transit, imposes this asymmetry among the covariant parameters i?s, -Rp and b. 

Our downward revision of the planetary radius translates into an increased value for the 
mean density, 0.67 ± 0.07 gem -3 . This value is 45-56% that of Jupiter. This is comparable 
to, but slightly less than, the mean densities of TrES-1 (0.84 gem -3 ; Sozzetti et al. 2004) 
and HD 189733b (0.93 gem" 3 ; Bakos et al. 2006). For XO-lb's estimated equilibrium tem- 
perature T eq = 1100 K (assuming Bond albedo A B = 0.4 and our derived value of the stellar 
radius R$ = 0.928 R Q ) and its mass Mp = 0.9 Mj, the models of Bodenheimer et al. (2003) 
predict planetary radii of Rp = 1.04 Rj and 1.11 Rj for models with and without a 20 M e 
core, respectively. Our estimate of the radius of XO-lb is 2 a larger than its predicted value, 
even for a planet without a core. HD 189733b's measured radius Rp = 1.154 ± 0.032 Rj 
(Bakos et al. 2006) is also larger than its theoretical value (Rp = 1.03 Rj with a core, 
Rp = 1.11 Rj without a core), given its equilibrium temperature T eq = 1050 K and mass 
M P = 0.82 ± 0.03 (Bouchy et al. 2005). In contrast, Laughlin et al. (2005) showed that 
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TrES-l's measured radius R P = 1.08 ± 0.05 Rj (Laughlin et al. 2005) is consistent with its 
theoretically predicted value (R P = 1.05 Rj with a core, R P = 1.09 Rj without a core). 

The measured radii of both XO-lb and HD 189733 are consistent with the predictions of 
Bodenheimer et al. (2003) if "kinetic heating" is included. In these models ~ 2% of the stellar 
insolation is deposited at depth, following the work of Guillot & Showman (2002). The age 
of XO-1 is also uncertain; the planet would be larger if the system were younger (Burrows 
et al. 2003). Whether or not the radius of XO-lb requires an additional energy source, as is 
the case for HD 209458b (0.36 gem -3 ; Knutson et al. 2006), is an important topic for future 
theoretical work. The kinetic heating model, proposed to explain the apparent inflation of 
HD 209458b, would naturally predict that many other "hot Jupiters" should be inflated. 
Other explanations, such as ongoing tidal circularization due to an eccentricity exchange 
with a third body (Bodenheimer et al. 2001), or the trapping in a Cassini state with nonzero 
obliquity (Winn & Holman 2005), would seemingly have difficulty accounting for a large 
population of inflated objects. 

The accuracy of our transit times ranges from 0.2 minutes for the FLWO z band observa- 
tions to 2.5 minutes for the transit observed solely by TopHAT. Fig. 5 shows the differences 
between the observed and predicted times of mid-transit, as a function of transit epoch. 
The predicted times assume the average orbital period determined by McCullough et al. 
(2006) and a reference time based on our observations. So far, all the times are marginally 
consistent with a constant period. These observations provide accurate anchors for future 
searches for transit time variations. 
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Fig. 5. — The timing residuals for the 4 observed transits, according to the ephemeris of 
McCullough et al. (2006). (see Eq. 1). The first point corresponds to the T c of McCullough 
et al. (2006). The points lie on a horizontal line, and therefore the data are marginally 
consistent with a constant period. 
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Table I. Photometry of XO-1 



Telescope 


Filter 


HJD 


Relative flux 


Uncertainty 


FLW048 


z 


2453883.70727 


1.00125 


0.00147 


Palomar60 


R 


2453875.77023 


1.00144 


0.00215 


TopHAT 


I 


2453875.64414 


1.00889 


0.00386 



Note. - The time stamps represent the Heliocentric Julian Date 
at the time of mid-exposure. The uncertainty estimates are based on 
the procedures described in § 2. We intend for this Table to appear 
in entirety in the electronic version of the journal. A portion is shown 
here to illustrate its format. The data are also available in digital 
form from the authors upon request. 
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Table 2. System Parameters of XO-1 



Parameter 


Median value 


68% conf. limits 






lower 


upper 


R s [R Q ] 


0.928 


— u.uoy 


+0.015 






n ni Qb 
— U.Ulo 


+U.Ulo 






— 0.032 


+0.034 


Rp [Rjup] 


1.184 


— 0.014 


+0.025 






— 0.018 


1 n nocb 
+0.028 






n ri/i oc 
— U.U4z 


+U.U4/ 


Rp/Rs 


0.13102 


— U.UUUD4 


+U.UUUD4 


b 


0.14 


-0.10 


+0.09 


i [deg] 


89.31 


-0.53 


+0.46 


tlV - h [hr] 


2.992 


-0.015 


+0.013 


tn — h [min] 


21.18 


-0.47 


+0.81 


m(z) 


0.128 


-0.071 


+0.061 


u 2 (z) 


0.60 


-0.14 


+0.12 


2ui(z) + u 2 (z) 


0.858 


-0.055 


+0.049 


ui(z) - 2u 2 {z) 


-1.07 


-0.34 


+0.30 


T c (17) [HJD] 


2453875.92305 


-0.00036 


+0.00032 


T c (18) [HJD] 


2453879.8640 


-0.0011 


+0.0010 


T c (19) [HJD] 


2453883.80565 


-0.00019 


+0.00017 


T c (20) [HJD] 


2453887.74679 


-0.00016 


+0.00014 



Note. — The parameter values in Column 2 are the median val- 
ues of the distributions shown in Fig. 3. The confidence limits in 
Columns 3 and 4 are based on the MCMC analysis. 

a These uncertainties ignore the uncertainty in stellar mass. (For 
parameters with no designation, the uncertainty in the stellar mass 
is irrelevant.) 

b These uncertainties include the 0.03 M uncertainty in the stellar 
mass reported by McCullough et al. (2006), propagated according to 

R s oc (Mg/M©) 1 / 3 and R P oc (Mg/M©) 1 / 3 . 

c These uncertainties include a stellar mass uncertainty of 0.10 M . 



